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Millimeter Imaging of the /3 Pictoris Debris Disk: 
Evidence for a Planetesimal Belt 

David J. Wilner^, Sean M. Andrews^, A. Meredith Hughes^'^ 

ABSTRACT 

We present observations at 1.3 millimeters wavelength of the (3 Pictoris debris disk 
with beam size 4'.'3 x 2'.'6 (83 x 50 AU) from the Submillimeter Array. The emission 
shows two peaks separated by ~ 7" along the disk plane, which we interpret as a highly 
inclined dust ring or belt. A simple model constrains the belt center to 94 ±8 AU, close 
to the prominent break in slope of the optical scattered light. We identify this region 
as the location as the main reservoir of dust producing planetesimals in the disk. 



Subject headings: circumstellar matter ■ 
stars: individual (/3 Pictoris) 



■ planetary systems: planet-disk interactions- 



Introduction 



^The discovery of excess far-infrared emis sion from the n earby (19.44 it 0.05 pc; Ivan Leeuwen 

2007 ) A6V-type main-sequence star Pictoris ( Au mann^^l98d) together with optical imaging of scat- 
tered light from circumstellar dust (jSmith fc Terril6,1984 ) established the "debris disk" paradigm 
where dust g rains orbitirig the star originate from an eroding reservo ir of larger bodies (see re- 



views by e.g. lArtvmowicz 



1997 



Backman &: Paresce 



surroundin g this young star (12^ ^ Myr; 



1993; 



Zuckerman et al 



2001 



WvattI I2OO8I ) . The nearly edge-on disk 



. „ ^ is relatively luminous {Fdisk/F* 

2.5 X 10~^: iLagrange et al.ll200d ) and has been studied in great detail with a pa.noply of observa- 



200c ; 



tional techniques. High resolution ima ges in the optical (lKa,las fc Jewitt ' 1995; Hea p et al 
Golimowski et al.ll2006l ') . near-infrared (IMquillet et al.lll997al : iTamura et al...2006. : .Boccaletti et al 



2009) and mid-infrared dwahhai et al.ll2003 : 1 Weinberger et al.i2003 : Okamoto et"aLil2004 - ll^lesco et al 



2OO5I ) all show a wealth of structure, including density concentrations, an inner cavity, and asymme- 
tries such as warps. These features, including a secondary disk of scattered light inclined by about 5° 
(jAhmic et al.ll2009l ). have been variously ascribed to the gravitational influence of a giant planet or 
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plane ts (e.g. iMouillet et al.lll997bl : lAugereau et al.ll200ll : iFreistetter et al.1 120071 : iKennedv Wvatt 
20ld ). Indeed, a plan etary mass companion at a p rojected distance of 8 AU from the star now has 
been directly imaged (jLagrange et al.ll2009l . l20ld ). 



The emerging view of debris disks like (3 Pictoris postulates the pr esence of a planetesimal 
belt that produces dust with a range of sizes through collisional cascades (jStrubbe fc: Chiangll200a : 
WvattI I2OO8I : iKuchner &: Starkll2010l ). The stirring of the planetesimals may be due to the gravity 
of ~ 1000 km-sized objects forme d within the belt (IKenvon &: Bromlevll2004l ). or to the presence of 
planets located closer to the star (jMoro-Martin et al.ll2007l ). In either case, the dynamical effects of 
stellar radiation create a distribution of grain sizes that depends on distance from the star, e.g. the 
blow-out of the smallest "/3-meteoroid" grains. An important consequence is that images of debris 
disks at different wave lengths are dominated by different grain sizes and can appear remarkably 
different (IWyatti l2006l ) . Observations at millimeter wavelengths are most sensitive to large grains 
that are minimally affected by radiative forces and thus have the potential to trace best the location 
of the dust producing parent planetesimals. The debris disk around Vega, for exa mple, shows a 



clunipy ring confined to radii < 200 AU at wavelengths of 350 /um and longward ([Holland et al 



1998 



Wilner et al.l 120021 : iMarsh et al.l |2006| ) , while it appears smooth and featureless and extends 



to radii ^ 800 AU in mi d-infrared light that arises predominantly from small grains expelled by 
radiation (jSu et al.l l2005l ). The debris disk around HR 8799, an A-type star that harbors t hree 
directly imaged planets, shows similar morphological changes with wavelength (ISu et al.ll2009l ). 



For the /3 Pictoris debris disk, the angular resolutions of (sub-)millimeter images have been 
too coarse to reveal much structure. Images from several different telescopes generally show dust 
emission extended along a po sition angle of 30° , consistent with the optical disk: JCMT/SCUBA 
at 850 Mm with a If beam (jHolland et al.lll998l ). APEX/LABOCA at 870 um. with an 18" beam 
dNilsson et al.lbood '). SEST/SIMBA at 1200 /xm with a 24" beam dLiseau et al.ll2003l). and Her- 



schel /SPIRE at 250, 350 and 500 /im with 18, 25, and 37" beams, respectively (jVandenbussche et al 
2OIOI ) . A separate peak or blob of dust emission is also found ^ 30" to the sou thwest of the star, but 



the re lationship of this peak to the main disk is unclear; iDent et al.l (|2000l ) and lVandenbussche et al 



(|2010l ) have suggested it is a background galaxy with a coincidental alignment with the disk plane. 



Millimeter interferometry offers a way to obtain higher angular resolution and more information 
on the largest detectable grain populations within the debris disk. In this Letter, we present imaging 
observations of (3 Pictoris at 1.3 millimeters wavelength from the Submillimeter Array (SMAfl that 



reveal a belt of emission around the star centered near a radius of 
reservoir of planetesimals. 



95 AU that likely marks a 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academ- 
ica Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academica 
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Observations 



We used the eight element SMA (|Ho et al.ll2004l ) on Mauna Kea, Hawaii to observe /3 Pictoris 



in the compact-north configuration (baseUnes 6-97 m) on 2010 August 9, and in the extended 
configuration (baseUnes 12-178 m) on 2010 September 1. The phase center was chosen to be 
a = 5^47'"17!09, 5 = -51°03'59'.'5 (J2000), about 0'.'9 from the stehar position at the epoch of the 
observations. The (3 Pictoris system is a chahenging target for the SMA as it never rises above 
20° elevation. Nonetheless, usable data were obtained in both configurations over the hour angle 
range ±1.7. The weather conditions were good on both days, with 225 GHz atmospheric opacities 
0.07-0.09 and stable atmospheric phase. The correlator was configured to provide the maximum 
4 GHz of bandwidth in each of two sidebands centered ±6 GHz from a central LO frequency of 
235.6 GHz (wavelength of 1.3 millimeters), with a uniform spectral resolution of 0.8125 MHz. At 
this frequency, the primary beam size is ~ 54" (FWHM). Observations of the strong source 3C454.3 
obtained at the start of each track were used to calibrate the passband response. Observations of 
the quasars J0538-440 and J0522-364 were interleaved with j3 Pictoris in order to calibrate time 
dependent gain variations. The astrometric uncertainty is < 0'.'3. The absolute flux scale was 
set with reference to observations of the standard calibrator Callisto in each track and should be 
accurate to better than 15%. The calibration procedure was performed using the IDL based MIR 
software. Subsequent imaging and deconvolution were done within the MIRIAD package. 



3. Results and Analysis 
3.1. 1.3 Millimeter Emission 

Figure [1] shows a contour image of the 1.3 millimeter emissio n overlaying a Hub ble Space Tele- 



scope/STIS coronographic image of optical scattered light from iHeap et al.l ((20001). The 1.3 mil- 
limeter image was made using natural weighting and a modest taper in the east-west direction 
to avoid extreme ellipticity of the synthesized beam, which is 4'.'3 x 2'.'6 (83 x 50 AU) oriented 
nearly north-south (position angle 2°). The maximum sidelobes of the dirty beam obtained with 
this weighting scheme are located about 13" to the east and west, with amplitude 15% of the 
central peak. The rms noise in this image is 0.6 mJy beam~^. The star symbol is plotted offset 
by (0"35, 0"70) from the phase center, within the uncertainties of the stellar position corrected for 
proper motion. The 1.3 millimeter emission shows two peaks at positions symmetrically offset from 
the stellar position by ^ 3"5 to the northeast and southwest, respectively. This basic morphology 
suggests a highly inclined ring or belt, where the peaks are due to limb brightening at the ansae 
(where the column density is highest). While the southwest peak appears slightly brighter, the 
difference lies within the noise and cannot be considered significant. 



Because the SMA observations are not sensitive enough to detect the stellar photosphere, 
the alignment of the images from the SMA and Hubble Space Telescope is limited by the absolute 
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Fig. 1. — SMA image of the 1.3 mill imeter continuum emission from (3 Pictoris overlayed on an 
image of optical scattered light from iHeap et al.l (|200d ). The contour levels are —2,2,4,6,... x 
0.6 mJy (the rms noise level). Negative contours are dotted. The ellipse in the lower left corner 
represents the 4'.'3 x 2'.'6 (FWHM) synthesized beam size. The star symbol indicates the location 
of the stellar photosphere. 



astrometry. Even taking account of this uncertainty, it seems that the two mil limeter emissioii peak s 
do not align perfectly along the 30° position angle of the primary optical disk (jKalas &: Jewittlll995l ). 
Instead, examination of Figure [1] suggests that the peaks align more closely with the 34° position 
angle of the scattered light secondary disk described bv iGolimowski et al.l (|2006l ). Observations 
with better resolution and sensitivity are needed to confirm this suggestion; the non-circular beam 
makes it difficult to assess small differences in orientation, and the millimeter emission structure 
itself may prove to be warped or more complex. 

The 1.3 millimeter flux in the detected structure is 13 it 1.4 mJy, estimated by integrating over 
the emission in the image. This value is on ly approximately h alf of the 24.3 it 3.0 mJy measured 
in the 24" SEST beam at 1.2 millimeters (jLiseau et al.l l2003l ) . a discrepancy significantly larger 
than expected from the mutual absolute calibration uncertainties and the spectral slope. The 
difference suggests the presence of an additional, extended 1.3 millimeter emission component, 
missed in these observations by the spatial filtering properties of the interferometer. Given the 
shortest SMA ba selines, the peak brig htness is diminished already by 50% for a 20" (FWHM) 
Gaussian source (jWilner &: Welchlll994l ). a size scale smaller than the SEST beam. Judging from 
the partially resolved images from far-infrared and submillimeter filled aperture telescopes, this 
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missing component is likely elongated along the disk, which extends beyond the SMA field of view. 

The separate dust peak to the southwest is detected at ~ 5a, offset by (-21'.'4±0'.'4, -22'/4±0'.'6) 
from the center of the field (not shown). The corresponding absolute position is a = 5'^47™14';82, 
6 = — 51°04'2l'.'9 (J2000). Because the primary beam correction is large and uncertain at this 
location beyond the half power point, it is difficult to provide an accurate estimate of the flux. 
The position is well determined, however, and shows that this peak does not lie along an extension 
of the optical disk. This supports previous suggestions that this feature is a background source, 
presumably a dusty galaxy. 

3.2. Belt Location and Width 

We constrain the basic properties of the 1.3 millimeter emission with a simple model that 
assumes the structure is characterized by a flat, axisymmetric belt of emission. We take the radial 
profile of the emission to be r~^'^, which is physically motivated by optically thin emission for 
constant surface density and a temperature profile of r~^'^ due to stellar irradiation. We fix the 
inclination and orientation of the belt on the sky to 87° and 33°, respectively; due to the limited 
resolution, small variations in these geometric parameters do not have significant effects on the 
results. This simple model has three parameters: the belt center R, belt width AR, and flux F. 
We calculate a grid of models over the parameter ranges 60 < i? < 130 AU and 2 < AR < 110 AU in 
steps of 2 AU, and 10 < F < 17 mJy in steps of 0.5 mJy, and calculate values for each model using 
all of the SMA visibilities. The right panel of Figure [2] shows the resulting surface (marginalized 
over the parameter F). The cross marks the best fit at = 94 ± 8 AU, AR = 34^32 AU, 
F = 15 lb 2 mJy; the uncertainties represent the formal la errors. The data strongly constrain the 
belt center location and allow for widths up to sizes comparable to the resolution of the observations. 
The three left panels of Figure [2] show the 1.3 millimeter image from Figure [1] together with images 
of the best-fit model and the residuals, all made in the same way. The model reproduces the main 
features of the data, and the residuals are consistent with noise. If a steeper radial emissivity were 
assumed, then the outer edge of the emission could extend further. However, a belt with a width 
that encroaches much closer to the star than the best fit may be difficult to reconcile with the 
mid-infrared emission from the system. A proper model that considers the constraints of the full 
spectral energy distribution requires many more assumptions than made here, in particular about 
the grain composition, grain size distribution, collisional behaviors and dynamics. 

4. Discussion 

The new millimeter observations improve substantially on previous single dish images and 
start to resolve fine structure in the /3 Pictoris disk. Since large grains cannot travel far from their 
place of origin due to short grain-grain collisional timescales and negligible radiation effects, the 
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Fig. 2. — Left panels: SMA image of the 1.3 millimeter emission from f3 Pictoris together with the 
image of the best-fit axisymmetric belt model and the residuals. The contour levels and beam size 
are the same as in Figured! The dashed line indicates a position angle of 34°. Right panel: The 
surface for the belt center and width model parameters, with contours at 1, 2, 3cj. The cross marks 
the best-fit model. 



emission at this long wavelength should trace the dust producing planetesimals. Inspection and 
analysis of the resolved millimeter emission suggest a highly inclined ring or belt centered at a 
radius within or near ~ 95 AU. Aside from the nearly edge-on viewing geometry, the millimeter 
morphology is strikingly simil ar to other well-st udied A-type stars wit h substantial circurn stellar 
dust, in particular Vega (e.g. iMarsh et al.l l2006l ) and Fomalhaut (e.g. [Holland et al.ll2003l ). The 



region interior to t his belt in Pictor i s is clearly not em pty, as evidenced by mid-infrared imaging 



and spectroscopy (jTelesco et al 



2005 



Chen et al.l 120071 . e.g.), but it must be relatively deficient in 



dust mass or millimeter-sized grains or both. 

The location of the millimeter emission belt corresponds closely to a prom inent break in the 
slope of the optical scattered hght, as well as a change in the optical color gradient (jGolimowski et al. 
20061 ). These properties are plausibly explained in a model with dust producing planetesimals lo- 
cated just interior the break , with stellar radiation (and possibly also a st ellar wind) creating a 



radial gradient in grain size (lAugereau et al. 



2001 



Strubbe fc Chiang) l2006l ). In this scenario, the 



extended halo of emission along the disk plane would be dominated by a population of small grains 
blown out onto highly elliptical or hyperbolic orbits, possibly with temperatures above the local 
blackbody values, that cover a large area on the sky and give rise to the fraction of millimeter 
emission missed by the interferometer. 

One implication of the multi-component emission structure is t hat the far- infrared to millim eter 
spectral index of 2.34 it 0.07 indicated by the integrated spectrum (jVandenbussche et al.ll2010l ) may 
not be representative of any of the individual components. Depending on the details, it is possible, 
for example, that the belt component could show a steeper spectral index that would be closer 
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to expectations for a steady-state collisional cascade, without resorting t o unusual fragmentation 



prescriptions or wavy grain size distributions (jXhebault &: Augereaull2007l ) . Until multi- wavelength 



observations are available in this regime that clearly resolve the relevant structures, it will be 
problematic to use the integrated spectral index to make conclusive inferences about the grain 
properties and size distributions. 

The 1.3 millimeter image of the dust belt around /3 Pictoris sets the stage for much improved 
future observations with the Atacama Large Millimeter/Submillimeter Array (ALMA), now under 
construction in Chile (and much better placed for aperture synthesis observations of this southern 
source). More detailed millimeter images have the potential to determine, e.g. if the dust belt 
center is offset from the star, or if the emission exhibits pericenter glow or other asymmetries that 
could point to dynamical perturbations from additional planets in this remarkable system. 
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